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ABSTRACT 



We calculate rapidity distribution of photons due to Compton and annihilation processes 
from Quark Gluon Plasma (QGP) with pre-equilibrium momentum-space anisotropy. We 
also include contributions from hadronic matter with late stage transverse expansion. A phe- 
nomeno logical model has been used for the time evolution of hard momentum scale Phard( r ) 
and anisotropy parameter £(t). As a result of pre-equilibrium momentum-space anisotropy, 
we find significant modification of photons rapidity distribution. For example, with fixed 
initial condition (FIC) free- streaming (5 = 2) interpolating model we observe significant 
enhancement of photon rapidity distribution at fixed Pt, where as for FIC collisionally- 
broadened (5 = 2/3) interpolating model the yield increases till y ~ 1. Beyond that sup- 
pression is observed. With fixed final multiplicity (FFM) free- streaming interpolating model 
we predict enhancement of photon yield which is less than the case of FIC. Suppression is 
always observed for FFM collisionally-broadened interpolating model. 



1 Introduction 

Relativistic heavy ion colliders (RHIC), at Brookhaven National Laboratory and upcoming 
Large Hadron Collider (LHC) at CERN, are designed to produce and study strongly inter- 
acting matter at high temperature and/or density. Experiments at the RHIC have already 
demonstrated that high px hadrons in central A + A collisions are significantly suppressed in 
comparison with that in binary scaled p + p collisions [1]. This observation has been referred 
to as jet-quenching which clearly indicates towards the formation of Quark Gluon Plasma 
(QGP) 3 at the RHIC experiments. Another most important task is to characterize different 
properties of this new state of matter, such as isotropization/thermalization. 

The most difficult problem lies in the determination of isotropization and thermalization 
time scales (r iso and r therm ) 4 . Studies on elliptic flow (upto about pt ~ 1.5 — 2 GeV) using 
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3 Theoretically QGP is expected to be formed when the temperature of nuclear matter is raised above its 
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critical value, T c ~ 170 MeV, or equivalently the energy density of nuclear matter is raised above 1 GeV/ fm 
4 From now on we will concentrate on the most simplest possibility that is both the time-scale are the 

Same, Ttherm = T"iso- 



ideal hydrodynamics indicate that the matter produced in such collisions becomes isotropic 
with T iso ~ 0.6 fm/c [2, 3, 4]. On the contrary, perturbative estimates yield much slower 
thermalization of QGP [5]. However, recent hydrodynamical studies [6] have shown that 
due to the poor knowledge of the initial conditions, there is a sizable amount of uncertainty 
in the estimation of thermalization or isotropization time. The other uncertain parameters 
are the transition temperature T c , the spatial profile, and the effects of flow. Thus it is 
necessary to find suitable probes which are sensitive to these parameters. Electromagnetic 
probes have long been considered to be one of the most promising tools to characterize the 
initial state of the collisions [7, 8]. Because of the very nature of their interactions with 
the constituents of the system they tend to leave the system without much change of their 
energy and momentum. Photons (dilepton as well) can be one of such observables. 

But, photons can carry information about the plasma initial conditions [9, 10, 11] only if 
the observed flow effects from the late stages of the collisions can be understood and modeled 
properly. The observation of pronounced transverse flow in the photon transverse momentum 
distribution has been taken into account in model calculations of photon p T distribution at 
various beam energies [9, 10, 12, 13, 14]. It is found that because of the transverse kick 
the low energy photons populate the intermediate regime and consequently, the contribution 
from hadronic matter becomes comparable with that from the hadronic matter destroying 
the window where the contribution from QGP was supposed to dominate [14]. Apart from 
transverse flow effects, the investigation of longitudinal evolution using HBT correlation 
measurements has been done in Ref. [15]. It is shown that the decrease of -R S ide with pr 
(> 2.5 GeV) provides a good indication whether transverse flow is significant or not. But, 
so far as the pre-equilibrium emission is concerned this effect is not important. However, in 
the late stage, transverse expansion becomes important and we include this while estimating 
the photon yield from the hadronic matter. 

Photon (dilepton) production from relativistic heavy-ion collisions has been extensively 
studied in Ref. [16, 17, 18, 19]. All these works are based on the assumption of rapid 
thermalization of plasma with r t h C rm = where r; is the time scale of plasma formation. 
However, due to rapid longitudinal expansion of the plasma at early time this assumption 
seems to be very drastic because it ignores the momentum space anisotropy developed along 
the beam axis. 

The phenomenological consequences of early stage pre-equilibrium momentum space 
anisotropy of QGP have been studied in Ref. [20, 21] in the context of dileptons and in 
Ref. [22, 23] in the context of photons. In Ref. [22] the effects of time-dependent momentum- 
space anisotropy of QGP on the medium photon production are discussed. It is shown that 
the introduction of early time momentum-space anisotropy can enhance the photon produc- 
tion yield significantly. Also the present authors calculate transverse momentum distribution 
of direct photons from various sources by taking into account the initial state momentum 
anisotropy of QGP and the late stage tranverse flow effects [23]. The total photon yield is 
compared with the recent measurement of photon transverse momentum distribution by the 
PHENIX Collaboration to extract the isotropization time [23]. It is found that the data can 
be reproduced with r iso in the range 0.5 — 1.5 fm/c. All these works show that the intro- 
duction of pre-equilibrium momentum space anisotropy has significant effect on the medium 
dilepton as well as photon productions. In the present work, we will be investigating the 
rapidity dependence of thermal photon in the presence of pre-equilibrium momentum space 
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anisotropy. The importance of rapidity distributions of photons and dileptons produced in 
relativistic heavy ion collisions has been previously realized in Ref. [24, 25, 26]. As for exam- 
ple, fluctuations in the rapidity distributions can signal a phase transition, a supercooling of 
the fluid, or the presence of quark-matter bubbles, etc. Moreover, it was shown in Ref. [15] 
that photon rapidity density can be a good probe to distinguish between Landau-like [27] 
and Bjorken-like [28] dynamics. Therefore, photon rapidity density should carry the signa- 
tures of pre-equilibrium momentum space anisotropy. The rapidity distribution of thermal 
photons produced in Pb + Pb collisions at CERN SPS is also demonstrated in Ref. [29] using 
a three-fluid hydrodynamical model. It is argued that rapidity dependence of photon spectra 
can provide cleaner insight about the rapidity dependence of the initial conditions, e. g. the 
temperature / time. 

In absence of any precise knowledge about the dynamics at early time of the collision, 
one can introduce phenomenological models to describe the evolution of the pre-equilibrium 
phase. In this work, we will use one such model, proposed in Ref. [20], for the time depen- 
dence of the anisotropy parameter, £(r), and hard momentum scale, Phard(' r )- This model 
introduces four parameter to parameterize the ignorance of pre-equilibrium dynamics: the 
parton formation time (r^), the isotropization time (r; so ), which is the time when the system 
starts to undergo ideal hydrodynamical expansion and 7 sets the sharpness of the transition 
to hydrodynamical behavior. The fourth parameter 5 is introduced to characterize the na- 
ture of pre-equilibrium anisotropy i.e whether the pre-equilibrium phase is non-interacting or 
collisionally broadened. The phenomenological model in Ref. [20] assumes Bjorken's boost 
invariant expansion of the plasma. Therefore, this leads to a rapidity independent produc- 
tion of thermal photons. One can easily see that this can not be true for collisions involving 
nuclei having a finite energy. The experimental data are expected to be better described 
by rapidity dependent parton distribution functions. Therefore, in this work, we have used 
rapidity dependent quark and anti-quark distribution functions. The rapidity dependence 
of the distribution functions arises from the rapidity dependence of the initial temperature, 

The organization of the paper is as follows. In the next section (section 2) we shall discuss 
the mechanisms of photon production rate using an anisotropic phase space distribution along 
with the space-time evolution of the matter. Results are presented in section 3 and finally 
we conclude in section 4. 

2 Formalism 

2.1 Photon rate : Anisotropic QGP 

The lowest order mechanisms for photon emission from QGP are the Compton (q(q) g — > 
9(9)7) an d the annihilation (qq — > #7) processes. The rate of photon production from 
anisotropic plasma due to Compton and annihilation processes has been calculated in Ref. [30] . 
The soft contribution is calculated by evaluating the photon polarization tensor for an oblate 
momentum-space anisotropy of the system where the cut-off scale is fixed at k c ~ y^Phard- 
Here phard is a hard-momentum scale that appears in the distribution functions. 

The differential photon production rate for 1 + 2 — > 3 + 7 in an anisotropic medium is 
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given by [30]: 



J« Mf d°p x d°p 2 d°p 3 

E d^ ~ 2^fl 2^ 1 (27r)32E 2 (27r)32i?3(27r)3 /l(Pl ' Phard(r '^' e)/2(P2 ' Phard(r '^' e) 

X (2ir) 4 5( Pl + P2 - P3 - P )\M\ 2 [1 ± / 3 (p3,Phard(T,?7),0] (1) 

where, \Ai\ 2 represents the spin averaged matrix element squared for one of those processes 
which contributes in the photon rate and M is the degeneracy factor of the corresponding 
process. £ is a parameter controlling the strength of the anisotropy with £ > —1. f\, f 2 and 
/ 3 are the anisotropic distribution functions of the medium partons and will be discussed 
in the following. Here it is assumed that the infrared singularities can be shielded by the 
thermal masses for the participating partons. This is a good approximation at times short 
compared to the time scale when plasma instabilities start to play an important role. 

The anisotropic distribution function can be obtained [31] by squeezing or stretching 
an arbitrary isotropic distribution function along the preferred direction in the momentum 
space, 

/,(p,e,Phard) = /HvV + ^P-n^PhardM) (2) 

where n is the direction of anisotropy. It is important to notice that £ > corresponds to 
a contraction of the distribution function in the direction of anisotropy and —1 < £ < 
corresponds to a stretching in the direction of anisotropy. In the context of relativistic heavy 
ion collisions, one can identify the direction of anisotropy with the beam axis along which the 
system expands initially. The hard momentum scale Phard is directly related to the average 
momentum of the partons. In the case of an isotropic QGP, phard can be identified with the 
plasma temperature (T). 

2.2 Photon rate: Hadronic matter 

Photons are also produced from different hadronic reactions from hadronic matter either 
formed initially (no QGP scenario) or realized as a result of a phase transition (assumed to 
be first order in the present work) from QGP. Photons from hadronic reactions and decays 
cannot be calculated in a model-independent way. The hadronic matter produced in heavy 
ion collisions is usually considered to be a gas of the low lying mesons n, p, u, r\ and nucleons. 
Reactions between these as well as the decays of the p and u were considered to be the sources 
of thermal photons from hadronic matter [7, 18, 32]. 

We follow the calculations done in Ref. [33] where convenient parameterizations have 
been given for the reactions considered. These parameterizations will be used while doing 
the space-time evolution to calculate the photon yield from meson-meson reactions. 

2.3 Space time evolution 

The rate given in Eq. (1) is the static rate which has to be convoluted with the space- 
time history of the plasma to obtain phenomenologically predictable quantities, for example, 
dN '/ 'cPpTdy for a given p T or dN/d 2 p T dy for a given y. In our calculation, we assume an 
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isotropic plasma is formed at initial temperature Tj and initial time (proper) Tj. Subsequent 
rapid expansion of the matter along the longitudinal direction causes faster cooling in the 
beam direction than in the transverse direction. As a result the system becomes anisotropic 
and remains so till r = r iso . 

The exact dynamics at the early-stage of the heavy ion collision is almost unknown. Thus, 
a precise theoretical picture of the evolutions of Phard and £ is not possible. However, we can 
always introduce phenomenological models to parameterize the ignorance. In this work, we 
shall closely follow the work of Ref. [21] to evaluate the rapidity distribution of photons from 
the first few Fermi of the plasma evolution. Three scenarios of the space-time evolution (as 
described in Ref. [21]) are the following: (i) r iso = r,, the system evolves hydrodynamically 
so that £ = and Phard can be identified with the temperature (T) of the system (till date 
all the calculations have been performed in this scenario), (ii) r iso — > oo, the system never 
comes to equilibrium, (iii) r iso > Tj and r iso is finite, one should devise a time evolution 
model for £ and Phard which smoothly interpolates between pre-equilibrium anisotropy and 
hydrodynamics. We shall follow scenario (iii) (see Ref. [21] for details) in which case the 
time dependence of the anisotropy parameter £ is given by 



where the exponent 5 = 2 (2/3) corresponds to free- streaming (collisionally-broadened) pre- 
equilibrium momentum space anisotropy and 5 = corresponds to complete isotropization. 
As in Ref. [21], a transition width 7 _1 is introduced to take into account the smooth transition 
from non-zero value of 5 to 5 = at r = r iso . The time dependence of various quantities are, 
therefore, obtained in terms of a smeared step function [20]: 



For r << T iso (>> T iso ) we have A = 0(1) which corresponds to free- streaming (hydrodynam- 
ics). With this, the time dependence of relevant quantities are as follows [21, 26]: 



Ti 



(3) 




(4) 




Phard(r,r?) 



TM W V3 (r), 



(5) 



where, 




(6) 



and Ti is the initial temperature of the plasma. 
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Figure 1: (Color online) Photon transverse momentum (pt) distributions at RHIC energy 
with initial condition Tj = 0.446 GeV and r, = 0.147 fm/c. MM (BM) represents the 
contribution from meson- meson (baryon- meson) interactions (see [23] for details). 

To estimate the initial conditions we assume that i.e the longitudinal expansion ap- 
proximately follows the scaling law, v z = z/t, we can relate the initial density to the final 
multiplicity distribution by 

1 dN 

U i T i = p 2 , ( 7 ) 

7rRj_ drj 

and the initial temperature (Tj) is related with the initial density by the following relation, 
where #q is degeneracy factor. The multiplicity distribution can be parameterized as 



drj \ drj I „ V 2cr 2 



where (dN/dr]) Q is the total multiplicity at rj = and a = 3(5) for RHIC (LHC) energies. 
The initial temperature is therefore a function of r], i. e Tj = Ti{rj). 

So far, we have discussed the evolution during early stage only. For (r > r iso ) we propose 
that r iso onward the system is described by (1 + 2)d ideal hydrodynamics. As the system 
becomes isotropic at r = r iso , Phard(' r iso) and r iso can be identified as the initial conditions, i. 
e., initial temperature and initial time for the hydrodynamic evolution. The initial conditions 
for ideal hydrodynamics is given by, 

yhydro / \ 

- 1 % /''hard V ' iso/ 

rf ydm = r iso (10) 
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Figure 2: (Color online) The distribution of photons for FIC free- streaming (5 = 2) 
interpolating model at RHIC energy at (a) y — and (b) y = 4. 

Eventually the system undergoes a phase transition (assumed to be first order in the present 
work) at r = r/, where 7/ is determined by the condition phard(T = 7/) = T c , where T c ~ 170 
MeV. The phase transition ends at r H = r d Tj, where r d = gq/gn is the ratio of the degrees 
of freedom in the two (QGP phase and hadronic phase) phases. Therefore, the total photon 
yield, arising from present scenario is given by the following equation, 



dN 



d 2 pTdy 



d xE—- 

d 4 p 



+ 



d xE—- 

d 6 p 



(11) 

hydro 



where the first term denotes the contribution from the anisotropic QGP phase and the second 
term represents the contributions evaluated in ideal hydrodynamics scenario. The rapidity 
density is defined as, 

dN_ = j d 2 pT _dN_ ^ 
dy J dyd 2 px 

In order to numerically compute the integrals in Eqs. (11) and (12), one needs to know 
the time dependence of Phard and £ which has been discussed earlier. 

After introducing the space-time evolution of relevant quantities, we are now completely 
equipped for the integration of Eqs. (11) and (12). However, before going into the details of 
rapidity distributions of photons at the heavy ion collider experiments like RHIC and LHC, 
let us concentrate on few interesting features of the phenomenological model introduced for 
the evolution of Phard(T, v) an d £( r )- 

The simple model, introduced in section 2.3, which smoothly interpolates between an 
initially non-equilibrium plasma to an isotropic plasma, is based on the assumption that 
the initial conditions are held fixed. The smooth interpolation (keeping the initial condition 
fixed) between anisotropic and isotropic phases (described in section 2) results into a hard 
momentum scale (for r >> r iso ) which is by a factor [TZ((r iso j 't.^ — l)] 0,25 larger compared to 
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Figure 3: (Color online) Rapidity distribution (y) of photon at RHIC energy from quark 
matter (QM) and hadronic matter (HM) for (a) px = 2 GeV and (b) pt = 3 GeV. 

the momentum scale results from the hydrodynamic expansion of a system (with the same 
initial condition) from the beginning (see Eqs. (5) and (6)). As a consequence of this en- 
hancement of Phard(T, ??), the fixed initial condition interpolating models (both free- streaming 
and collisionally-broadened) will result in generation of particle number during the transition 
from non-equilibrium to equilibrium phase. Moreover, the entropy generation increases with 
the increasing value of r iso . Thus, the requirement of bounded entropy generation can be used 
to put some upper bound on the value of r iso for fixed initial condition interpolating mod- 
els. As for example, if we allow maximum 20% entropy generation at RHIC, the maximum 
possible value of ri so will be 1.2 (18) fm/c for fixed initial condition collisionally-broadened 
(free- streaming) interpolating model. 

Due to the phenomenological constraints on the entropy generation, one might not allow 
any entropy generation at all. In that case, one can redefine U{r) in Eq. (5) to ensure 
fixed final multiplicity in this model. Since we know the amount of enhancement (which is 
respondable for this entropy generation) of Phard, the redefinition of U{r) will be straight 
forward [21]: 

U{t) = U{t) [ft((r iso /r^ - 1)] _3/4 (n/r^) (13) 

It is important to mention that this redefinition corresponds to a lower initial temperature 
(Phard(Ti, v) < ^i) f° r r iso > T i- Larger value of isotropization time corresponds to lower 
initial temperature. 



2.4 Photon rapidity distribution for fixed 

The rapidity distribution of photons at fixed pt from a QGP or hadronic matter can be 
obtained by integrating Eq. (11). We use a Monte-Carlo computer code to numerically 
evaluate Eq. (11). For RHIC energies the formation time is taken as r< = 0.147 fm/c. 
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Figure 4: (Color online) Photon rapidity distribution at RHIC energy for FIC interpolating 
model, at p T = 2 GeV for (a) 5 = 2 and (b) 5 = 2/3. 



This corresponds to an initial temperature of Tj = 446 MeV at rj = [19]. For the LHC 
energies, we assume the formation time to be T{ = 0.073 fm/c and the initial temperature to 
be Ti = 897 MeV at r] = 0. 

As a consequence of pre-equilibrium momentum-space anisotropy, significant modification 
of the thermal photon rapidity distribution is expected. To quantify the effect of isotropiza- 
tion time on the rapidity distribution of the thermal photons, we define photon modification 
factor $(?/, T iso )| PT for fixed pt as the ratio of photon yields with and without pre-equilibrium 
momentum-space anisotropy, 

/ dN(y, r iso ) \ / dN(y,T iso = n) \ 



The modification factor $ is not measurable quantity. If there is an anisotropy in nature, the 
numerator can be measured (see Eq. 14). If there is no anisotropy in nature the denominator 
can be measured. But it is impossible to measure both simultaneously. However, $ is a useful 
quantity only for demonstrative purpose. 



3 Results 

Before going to the description of photon rapidity distribution, we validate the present model 
in the context of PHENIX photon data [35]. It has already been demonstrated in Ref. [23]. 
We shall repeat here for completeness. 

To show the presence of initial state momentum anisotropy, we plot the total photon yield 
assuming hydrodynamic evolution from the very begining as well as with finite Ti SO in Fig. 1. 
It is clearly seen that some amount of anisotropy is needed to reproduce the data. We note 
that the value of r iso needed to describe the data lies in the range 0.5 fm/c< r iso < 1.5 fm/c, 
independent of values of the transition temperatures [23]. It is to be noted that observables 
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Figure 5: (Color online) Photon rapidity distribution at RHIC energy for FIC interpolating 
model, at p T = 4 GeV for (a) 5 = 2 and (b) 5 = 2/3. 

like photon px distribution along with the ratio of photon to pion rapidity density will further 
strengthen the validity of the present model. 

3.1 Fixed initial condition (FIC) interpolating model 

Fixed initial condition (FIC) interpolating models always result into an enhanced value of 
hard momentum scale as a consequence of pre-equilibrium anisotropy. This feature of this 
model has already been discussed briefly before (see Ref. [21] for details). As a consequence of 
this enhancement of Phard> pre-equilibrium anisotropy increases the density of plasma partons 
with small rapidities. However, for higher rapidities, this enhancement is complemented by 
the suppression arising from the non-zero value of anisotropy parameter £. The anisotropic 
parton distribution function in Eq. (2) clearly suggests that for the positive values of £, 
parton density decreases if we decrease the angle between the momentum of partons and 
the direction of anisotropy. From the very beginning, we have assumed that pre-equilibrium 
momentum-space anisotropy (in the heavy ion collisions) results from the rapid longitudinal 
expansion of the system immediately after the collision. Thus, in the context of relativistic 
heavy ion collision, one can identify the direction of anisotropy as the beam axis. Moreover, 
due to the rapid longitudinal cooling (as a consequence of rapid longitudinal expansion), 
one always finds oblate anisotropic distributions i.e positive value of £. Therefore, intro- 
duction of pre-equilibrium anisotropy with fixed initial condition increases the density of 
plasma partons moving in the transverse direction [22] and at the same time decreases the 
density of plasma partons moving in the forward direction. This feature of fixed initial condi- 
tion pre-equilibrium momentum-space anisotropy should be reflected in the photon rapidity 
distribution which we will discuss in the following. 

As a consequence of this enhancement of partons moving in the transverse direction and 
suppression of partons moving in forward direction, we expect an enhancement of the yield 
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Figure 6: (Color online) Modification factor for FIC interpolating model at RHIC energy, 
(a) p T = 2 GeV and (b) p T = 4 GeV. 

in the transverse direction (y = 0) 5 and suppression in the longitudinal direction (y ^ 0). 

To see the relative importance of the contributions from quark matter (QM) and hadronic 
matter (HM) we evaluate the pt distribution of photons for fixed rapidity from the evoluting 
fireball. Because the effect of tranverse flow is pronounced in the late stages of the collisions 
we neglect this effect in the early stage and include it during the late stage. For r > r iso , 
the system is described by ideal relativistic hydrodynamics in (1 + 2)d [36] with longitudinal 
boost invariance [28] and cylindrical symmetry. Therefore, the total thermal photon yield 
for fixed y, arising from the present scenario (FIC interpolating model) is given by Eq. (11). 

To elucidate the effect of transverse expansion, we first show results in the frame work of 
FIC free- streaming (5 = 2) interpolating model. In Fig. 2 we have plotted the pt distribution 
of photons for two different isotropization time (t[ SO = 1 and 2 fm/c) at (a) y = and (b) 
y = 4. Both Figs. 2a and 2b clearly suggest that for 1 < px < 4 GeV, the photons 
from the quark matter (QM) dominates over the hadronic matter (HM). As we increase the 
isotropization time (ri SO ) from 1 to 2 fm/c (see Fig. 2), both the contributions from quark 
matter and hadronic matter increases. However, quark matter contributions still dominate 
over hadronic matter contributions. Therefore, for the above mentioned pt range, the effects 
of pre-equilibrium momentum space anisotropy will be much more prominent irrespective 
of isotropization time (ri SO ). In support of this argument, we have plotted dN / cPpxdy as 
a function of rapidity (y) for fixed pt (pt — 2 GeV) in Fig. 3a for two different values of 
isotropization time (r iso = 1 and 2 fm/c respectively). Fig. 3a shows that the contributions 
from hadronic matter are always smaller than the contributions from the quark matter. 
However, for pt = 3 GeV, QM and HM contributions are of similar magnitude for y < 1 
(see Fig. 3b). Beyond that the HM contributions dominates. 

The feature described above is better understood by looking at Figs. 4 and 5 where the 
total contributions have been plotted. In Fig. 4, we have presented the rapidity distribution 

5 This feature was already established in our previous work [22]. 
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Figure 7: (Color online) Modification factor for FFM interpolating model at RHIC energy, 
(a) p T = 2 GeV,(b) p T = 4 GeV. 

of photons (diV/ d 2 pxdy) as a function of rapidity (y) with fixed transverse momentum = 
2 GeV (for three different values of isotropization time, Ti SO = Tj, 1 and 2 fm/c) in the 
framework of fixed initial condition for (a) free- streaming (5 = 2) and (b) collisionally- 
broadened (5 = 2/3) interpolating models. The rapidity distribution of photons with p? = 4 
GeV are presented in Fig. 5 for (a) free- streaming (5 = 2) and (b) collisionally- broadened 
(5 = 2/3) interpolating model. Figs. 4a and 5a show enhancements of photon yields (for 
the FIC free-streaming (5 = 2)) interpolating model) in the low rapidity region (0 < y < 5). 
Marginal suppressions in the higher rapidity region (y > 5) are observed. The situation is 
different for the FIC collisionally broadened (5 = 2/3) interpolating model (see Figs. 4b and 
5b) which shows slight enhancements of photon yields in the low rapidity region (0 < y < 1.5) 
and suppressions for the rest of the rapidity region (y > 1.5). This can be attributed to the 
fact that in the case of collisionally-broadened interpolating model we have included the 
possibility of momentum space broadening of the plasma partons due to interactions. As a 
consequence the hard momentum scale Phard (which is related to the average momentum in 
the partonic distribution functions) decreases with time, whereas for free- streaming model, 
the hard momentum scale remains unchanged (phard(' r ) = Ph&rd( T i) = T^, for r < r iso ) upto 

7~ 7"iso- 

These enhancement and suppression are more clearly revealed from Fig. 6 which shows 
the modification factors for (a) pt = 2 GeV and (b) px = 4 GeV. For a better comparison 
of free- streaming and collisionally-broadened pre-equilibrium phases, we have plotted the 
modification factors for 5 = 2 and 2/3 in the same graph. 

We always find enhancement (one order of magnitude) for 5 = 2. The trend of the graph 
shows that at larger rapidities (y > 5) the yield will be suppressed for the FIC free- streaming 
(5 = 2) interpolating model. However, we observe suppression for FIC collisionally-broadened 
interpolating model at rapidities y > 1.5. These findings are similar to that in Ref. [22] 
where for 5 = 2 with FIC considerable enhancement is observed at y = 0. But for 5 = 2/3, 
it is found that the enhancement is small. 
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Figure 8: (Color online) Modification factors, (a) FIC and (b) FFM interpolating models at 
LHC energy at pr = 4 GeV. 

3.2 Fixed final multiplicity (FFM) interpolating model 

The problems, regarding the entropy generation (appears in fixed initial condition model), 
could be eliminated by enforcing fixed final multiplicity. However, enforcing fixed final mul- 
tiplicity corresponds to a lower value of initial hard momentum scale (phard(^) < Ti{rfj). 
Moreover, in this case, the initial hard momentum scale will be isotropization time depen- 
dent i.e larger the value of t iso lower will be the initial hard momentum scale. The suppression 
of hard momentum scale corresponds to a suppression in the plasma parton density com- 
pared to the fixed initial condition interpolating model. Therefore, for fixed final multiplicity 
interpolating model, we predict more suppression in the photon yield compared to the FIC 
interpolating model. As a consequence of non-zero positive anisotropy parameter £, photons 
with higher rapidities will be more suppressed. 

These features (compared to fixed initial condition interpolating model) can be better 
understood by looking at the modification factors for fixed final multiplicity interpolating 
model displayed in Fig. 7 for (a) pt = 2 GeV and (b) pt = 4 GeV at RHIC energy. For 
5 = 2 (see Fig. 7a) enhancement is observed for the rapidity range upto y < 4 GeV. After 
that (y > 4), we observe the suppression of the photon yield. But for 5 = 2/3 (see the same 
Fig), we observe slight enhancement of the photon yield upto y < 1.5 and suppression for 
the rest of the rapidity region. For 5 = 2/3 at lower rapidities the rate is close to unity and 
as we go to higher rapidities the rate is reduced. It is to noted that although, in the case of 
FFM interpolating model, we have started with a lower value of initial hard momentum scale 
(phard(Ti) < T^Ti)), we obtain an enhancement in the low rapidity region (Fig. 7a) at px = 2 
GeV for both 5 = 2 and 2/3 respectively. For higher p T (4 GeV) the yield is suppressed at 
all rapidities irrespective of the values of 5 (Fig. 7b). However, the amount of suppression, 
for 5 = 2/3 at rapidities y > 2, is large compared to the free- streaming pre-equilibrium phase 
(5 = 2). For a more realistic case (5 = 2/3) the suppression is by a factor of 0.85 which is in 
accordance with Ref. [26] for the case of dilepton. 
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Figure 9: (Color online) Photon rapidity density at (a) RHIC and (b) LHC energies for FFM 
free-streaming (5 = 2) interpolating models. 



In Fig. 8, we have presented the photon modification factor for px = 4 GeV as a function 
of photon rapidity for ^snm — 5.5 TeV. Fig. 8a shows enhancement (suppression) for free- 
streaming (collisionally-broadened) pre-equilibrium phase for FIC. On the other hand, for 
FFM model (see Fig. 8b), we obtain suppression for collisionally-broadened pre-equilibrium 
phase. For free- streaming pre-equilibrium phase, we find small enhancement in the lower 
rapidity region and suppression for higher rapidities. 



3.3 Photon rapidity density 

There is one more phenomenologically interesting observable namely the photon rapidity 
density, dN/dy, which is defined in Eq. (12). Its shape at least is not expected to be affected 
by the transverse expansion of the system. Therefore, this quantity is phenomenologically 
very relevant. We have assumed two different values of isotropization time (t; so = 1,2 fm/c) 
to compute the rapidity density of photons from an anisotropic QGP. Since the detection of 
very soft photons is experimentally challenging, photons having their transverse momentum 
greater than 1 GeV is considered. 

The consequences of pre-equilibrium momentum-space anisotropy (with fixed initial con- 
dition and fixed final multiplicity) on the photon production have already been discussed in 
section 3.1 and 3.2 in details. 

We show the effects of tranverse expansion to the rapidity density in Fig. 9, where we 
have plotted the individual contributions of photon rapidity density from the quark matter 
(QM) and the hadronic matter (HM) in the frame work of free- streaming (5 = 2) FFM 
interpolating model at (a) RHIC and (b) LHC energies. As we vary t 1so from 1 to 2 fm/c 
(in both the Figs. 9a and 9b), the individual contributions from the quark matter and the 
hadronic matter increase. At the same time rapidity density of photons from the quark 
matter is higher than the contributions from the hadronic matter. This is true for the RHIC 
as well as LHC energies (see Fig. 9). This result is consistent with Fig. 3. 
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Let us now mention the main features of our observations. With FIC for both values 
of 5 we see moderate enhancement at not too large rapidities. Similar feature has been 
noted earlier [22]. This can be attributed to the fact that momentum anisotropy enhances 
the density of plasma partons in the tranverse direction. However, it is also to be noted 
that for larger rapidities the rate is suppressed for both the values of 5. This is due to the 
fact that the density of partons with higher longitudinal momentum (larger y) decreases. 
It is worthwhile to mention that the effects of pre-equilibrium momentum space anisotropy 
for 5 = 2 is more compared to that for 5 = 2/3. This is because 5 = 2/3 corresponds 
to close to isotropization. With FFM we observe marginal enhancement at low rapidities 
but considerable suppression at higher rapidities for both free- streaming and collisionally- 
broadened interpolating models. This suppression can be described in two ways. Due to 
rapid longitudinal expansion the distribution function becomes anisotropic. Photons with 
the larger values of longitudinal momentum are reduced compared with the photons with 
isotropic distribution function. Maximum amount of momentum-space anisotropy achieved 
in the early times will be the important cause of the suppression. The suppression will 
also depend on the time dependence of the anisotropy parameter £. Another source of 
rapidity dependence is given by hard momentum scale (phard) which depends on the initial 
temperature (Tj) and hence on rapidity. So we see that the hard momentum scale is directly 
related with the 77 even in the case of instantaneous thermalization and satisfies the relation 
in Eq. (5). 

It is important to mention that photon rapidity density significantly depends on the ex- 
pansion dynamics of the system. In absence of any theoretical knowledge, one can introduce 
different models for the expansion dynamics of the system like, Bjorken and Landau dynam- 
ics etc. In Ref. [15], it is shown that different expansion scenarios predicts different shape 
for photon rapidity density and thus, the shape of photon rapidity density can be used to 
distinguish different expansion scenarios. It is argued that the actual expansion scenario lies 
between Bjorken and Landau hydrodynamics by using various photonic observables [15]. 

In this article, we have only concentrated on the thermal photon rapidity (from both the 
quark matter and the hadronic matter) density as a probe of pre-equilibrium anisotropy. 
However, there are other non-thermal sources of photon like, hard photons, fragmentation 
photons, decay photons etc. Hard photons are produced in the initial hard scattering of col- 
liding nuclei and are insensitive to the later stage evolution of the QGP. Therefore, they do 
not carry any information about the pre-equilibrium anisotropy. Similarly, fragmentation and 
decay photons are also not sensitive to the pre-equilibrium phase. It is to be mentioned that 
apart from thermal photons from quark matter, thermal photons are also produced from hot 
hadronic matter where the late stage transverse expansion plays important role. In this work 
we include this effects of late stage transverse expansion. There is another important source 
of photons namely, the photons from jet-plasma interaction. The importance of jet-plasma 
interaction in the context of PHENIX [35] photon data was shown in Ref. [19]. For jet-plasma 
interaction, the jet-parton distribution functions do not depend on the pre-equilibrium phase. 
However, the plasma parton densities depend on the pre-equilibrium anisotropy. Therefore, 
the jet-plasma interaction is sensitive to the pre-equilibrium anisotropy. However, the sensi- 
tivity is smaller compared to the thermal photons. Because, for thermal photons, both the 
initial state partons are sensitive to the anisotropy. Moreover, the jet-photon contribution 
dominates over the thermal photon contribution only in the high region. Therefore, pho- 
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Figure 10: (Color online) thermal photon and jet-photon rapidity density at RHIC energy 
for FIC free- streaming (5 = 2) interpolating model. 

tons from jet-plasma interaction insignificantly contributes to the photon rapidity density in 
the pt range considered here. 

In Fig. 10, we have plotted the contribution from jet plasma interaction along with the 
anisotropic contributions for two values of isotropization times in the frame work of FIC free- 
streaming interpolating model. We have also presented the contribution from isotropic QGP 
(i.e. corresponding thermal photon rapidity density) in the same figure for comparison. 
It is seen that the contribution from jet-plasma interaction is well below the anisotropic 
contribution, leaving behind a window - where the effects of anisotropy could be seen. Fig. 
10 clearly shows that jet-photon contribution is not very significant in the context of thermal 
photon rapidity density. 



4 Conclusion 

To summarize, we have investigated the effects of the pre-equilibrium momentum space 
anisotropy of the QGP along with the effects of late stage tranverse expansion on the Pt{u) 
distribution at fixed y{pr) and rapidity density (dN/dy) of photons. To describe space-time 
evolution of hard momentum scale, Phard(T) and anisotropy parameter, £(r), phenomenolog- 
ical models have been used [21]. These phenomenological models assume the existence of 
an intermediate time scale called the isotropization time (r iso ). The first model is based on 
the assumption of fixed initial condition. However, enforcing fixed initial condition causes 
entropy generation. Therefore, we have also considered another model, which assumes the 
fixed final multiplicity. Both the possibilities of free- streaming and collisionally-broadened 
pre-equilibrium phase of the QGP are considered. The rapidity distribution of photons 
for different isotropization times in the frame work of these phenomenological models have 
been estimated. We observed that, for fixed initial condition, a free streaming interpolat- 
ing model can enhance the photon yield significantly for rapidities upto y ~ 4.5. However, 
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for collisionally-broadened pre-equilibrium phase with fixed initial condition, the enhance- 
ment of photon yield is upto y ~ 1.5. After that we observe the suppression of the photon 
yield for the entire rapidity region (y > 1.5). Since fixing the final multiplicity reduces the 
initial hard momentum scale or equivalently the initial energy density, we observe slight en- 
hancement in the low rapidity region and significant suppression (both for the free-streaming 
and collisionally-broadened interpolating models) for the rest of the rapidity region. This 
suppression can be explained as a consequence of the combined effect of the anisotropy in 
momentum-space achieved at early times due to expansion and the rapidity dependence of 
the hard momentum scale. For RHIC energies at pt = 2 GeV, QM contribution dominates 
over the HM contribution. However, for pt = 3 Gev, the later dominates over QM for y < 1. 
But as far as the total contribution is concerned, we always find significant modification 
(enhancement or suppression depending upon the initial conditions used) of the yield in 
presence of pre-equilibrium momentum space anisotropy. 

The other observables like heavy-quark transport [37], jet-medium- induced electromag- 
netic and gluonic radiation could be phenomenologically very useful in order to detect the 
consequences of pre-equilibrium momentum-space anisotropy. 

Acknowledgments: We would like to thank D. K. Srivastava and S. Sarkar for their 
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